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Abstract: We compute the screening length for quarkonium mesons moving through an 
anisotropic, strongly coupled = 4 super Yang-Mills plasma by means of its gravity 
dual. We present the results for arbitrary velocities and orientations of the mesons, as well 
as for arbitrary values of the anisotropy. The anisotropic screening length can be larger 
or smaller than the isotropic one, and this depends on whether the comparison is made 
at equal temperatures or at equal entropy densities. For generic motion we find that: (i) 
mesons dissociate above a certain critical value of the anisotropy, even at zero temperature; 
(ii) there is a limiting velocity for mesons in the plasma, even at zero temperature; (iii) in 
the ultra-relativistic limit the screening length scales as (1 — v'^Y with e = 1/2, in contrast 
with the isotropic result e = 1/4. 
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1 Introduction 



A remarkable conclusion from the experiments at the Relativistic Heavy Ion Collider 
(RHIC) [1, 2] and at the Large Hadron Collider (LHC) [3] is that the quark-gluon plasma 
(QGP) does not behave as a weakly coupled gas of quarks and gluons, but rather as a 
strongly coupled fluid [4, 5]. This renders perturbative methods inapplicable in general. 
The lattice formulation of Quantum Chromodynamics (QCD) is also of limited utility, since 
for example it is not well suited for studying real-time phenomena. This has provided a 
strong motivation for understanding the dynamics of strongly coupled non-Abelian plas- 
mas through the gauge/string duality [6-8] (see [9] for a recent review of applications to 
the QGP). 

For a period of time immediately after the collision, the system thus created is 
far from equilibrium. After a time Ti^^ > T^ut the system becomes locally isotropic and 
a standard hydrodynamic description becomes applicable. It has been proposed than an 
intrinsically anisotropic hydrodynamical description can be used to describe the system at 
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intermediate times r^nt < t < Ti^^ [10-18]. In this phase the plasma is assumed to have 
significantly unequal pressures in the longitudinal and transverse directions. The standard 
hydrodynamic description is a derivative expansion around equal pressures, and therefore 
it is not applicable in this regime. In contrast, the intrinsically anisotropic hydrodynamical 
description is a derivative expansion around an anisotropic state, and hence in this case the 
requirement that derivative corrections be small does not imply small pressure differences. 
In a real collision the degree of anisotropy will decrease with time, but for some purposes 
it is a good approximation to take it to be constant over an appropriate time scale. 

Motivated by these considerations, in this paper we will investigate the effect of an 
intrinsic anisotropy on the screening length between a quark-antiquark pair in the plasma. 
We will consider the case in which the pair is at rest in the plasma as well as the case in 
which it is moving through the plasma. For this purpose we will examine a string with 
both endpoints on the boundary of an asymptotically AdS spacetime [19, 20] that is dual 
to an anisotropic = 4 super Yang-Mills plasma. As we will review below, the plasma 
is held in anisotropic equilibrium by an external force. The gravity solution possesses an 
anisotropic horizon, it is completely regular on and outside the horizon, and it is solidly 
embedded in type IIB string theory. For these reasons it provides an ideal toy model in 
which questions about anisotropic effects at strong coupling can be addressed from first 
principles. 

To avoid any possible confusion, we clarify from the beginning that the quarks and 
antiquarks that we will consider are infinitely massive, i.e. the bound states that we will 
consider are the analogue of heavy quarkonium mesons in QCD. Thus, the reader should 
always have the word 'quarkonium' in mind despite the fact that we will often refer to these 
states simply as 'mesons', 'heavy mesons', 'quark-antiquark bound states', 'dipoles', etc. 
This is specially relevant in the ultra-relativistic limit of the screening length, to which 
we will pay particular attention since it can be determined analytically. We emphasize 
that our results correspond to sending the quark and antiquark masses to infinity first, 
and then sending t> — )• 1. In particular, this means that in any future attempt to connect 
our results to the phenomenology of the QGP, this connection can only be made to the 
phenomenology of heavy quarkonium moving through the plasma. 

The screening length for quarkonium mesons at rest in the anisotropic plasma of [19, 20] 
has been previously studied in [21, 22]. Our Sec. 4 has some overlap with these references 
and, wherever they overlap, our results agree with theirs. Other physical properties of the 
anisotropic plasma that have been calculated include its shear viscosity [23, 24], the drag 
force on a heavy quark [21, 25], the jet quenching parameter [21, 22, 26], and the energy 
lost by a rotating quark [27]. The phase diagram of the zero-coupling version of the model 
considered in [19, 20] has been studied in [28]. 
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Figure 1. Metric functions for a/T 
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4.4 (left) and a/T ~ 86 (right). 



2 Gravity solution 

The type IIB supergravity solution of [19, 20] in the string frame takes the form 

ds^ = ^ (-TB + dx^ + dy^ + ndz^ + —)+ L\¥d^l (2.1) 
X = az, ^ = Hu), (2.2) 

where x and cp are the axion and the dilaton, respectively, and {t, x, y, z) are the gauge 
theory coordinates. Since there is rotational invariance in the xy-directions, we will refer 
to these as the transverse directions, and to z as the longitudinal direction. J^, B and Ti 
are functions of the holographic radial coordinate u that were determined numerically in 
[19, 20]. Their form for two values of a/T is plotted in Fig. 1. The horizon lies at u = Un, 
where T = 0, and the boundary at w = 0, where T = B = H = 1 and (/> = 0. The 
metric near the boundary asymptotes to AdS^ x S^. Note that the axion is linear in the 
z-coordinate. The proportionality constant a has dimensions of mass and is a measure of 
the anisotropy. The axion profile is dual in the gauge theory to a position-dependent theta 
parameter of the form 6 cc z. This acts as an isotropy-breaking external source that forces 
the system into an anisotropic equilibrium state. 

If a = then the solution reduces to the isotropic black D3-brane solution dual to the 
isotropic = 4 theory at finite temperature. In this case 

1 

B = n = l, x = = O, F=l- — , u^ = ^ (2.3) 
and the entropy density takes the form 



s-... = -N^T\ (2.4) 

Fig. 2 (left) shows the entropy density per unit 3-volume in the xyz-directions of the 
anisotropic plasma as a function of the dimensionless ratio a/T, normalized to the en- 
tropy density of the isotropic plasma at the same temperature. At small a/T the entropy 
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Figure 2. (Left) Log-log plot of the entropy density per unit 3-volume in the xyz-directions as a 
function of a/T, with s^so defined as in eqn. (2.4). The dashed blue line is a straight line with slope 
1/3. (Right) Log- log plot of the entropy density per unit 2-area in the xy-directions on a constant-z 
slice as a function of a/T, with sf^ defined as in eqn. (2.6). 



density scales as in the isotropic case, whereas at large a/T it scales as [19, 20, 29] 

s = c^,,,N^a^/^T^/^ , [a/T > 1] (2.5) 

where c^nt is a constant that can be determined numerically. The transition between the 
two asymptotic behaviors of the entropy density takes place at a/T ~ 3.7. 

For later use in Fig. 2 (right) we show the entropy density per unit 2-area in the xy- 
directions on a constant-z slice, normalized to the isotropic result 

4° = ^N^T' . (2.6) 

Also for later use we list here the near-boundary behavior of the different functions that 
determine the solution (2.2): 

^ = 1 + ^a\2 + (^jr^ + log + 0(n6) , 

B = l- ^^^a'u^ + {b^ - ^a" log + ©(n^) , 

U = l + \a'u' - - - log tx) + 0{n^) . (2.7) 

The coefficients J-4 and B4 depend on a and T and are known analytically in the limits of 
low, and high temperature and numerically for intermediate regimes [20]. 

A feature of the solution (2.2) that played an important role in the analysis of [19, 20] 
is the presence of a conformal anomaly. Its origin lies in the fact that diffeomorphism 
invariance in the radial direction u gets broken in the process of renormalization of the on- 
shell supergravity action. In the gauge theory this means that scale invariance is broken 
by the renormalization process. One manifestation of the anomaly is the fact that, unlike 
the entropy density, other thermodynamic quantities do not depend solely on the ratio a/T 
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Figure 3. Ratio (2.8) as a function of a/T. The blue dots are the actual values of the ratio, and 
the red curve is the fit (2.9). 

but on o and T separately. Fortunately, this will not be the case for the screening length, 
as we will see below. 

To facilitate a (rough) comparison of the anisotropy in our system to that in other 
anisotropic plasmas it is useful to consider the ratio 



where E is the energy density and -Pl,Pl ai'e the transverse and longitudinal pressures, 
respectively. In addition to being dimensionless, this ratio has the virtue that it does not 
depend on a and T separately, but only on the combination a/T. For the isotropic A/" = 4 
super Yang-Mills plasma a = 1, whereas for < a/T < 20 the ratio is well approximated 
by the expression 



as shown in Fig. 3. 

At various points we will refer to the limit T = of the anisotropic plasma. The zero- 
temperature version of the solution (2.2) was found in [29]. In this case the string- frame 
metric exhibits a naked curvature singularity deep in the infra-red, and the Einstein-frame 
metric exhibits infinite tidal forces [30, 31]. However, we emphasize that, for any finite 
temperature, the singularity is hidden behind the horizon and the solution is completely 
regular on and outside the horizon, exhibiting no pathologies of any type. Thus we will 
think of the T = results as those obtained by taking the limit T — t- of the finite- 
temperature results. Moreover, regulating the infra-red geometry in this or any other way 
is actually unnecessary for most of the physics of quarkonium dissociation. The reason is 
that, as we will see, in the limit in which a/T becomes large the penetration depth into 
the AdS bulk of the string that is dual to the quarkonium meson becomes very small. As 
a result, the dissociation is entirely controlled by the metric near the boundary, which is 
insensitive to the infra-red behavior described above. 



a = 



AE + P^- Pi 
3Ts 
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(2.9) 
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Figure 4. Energy difference, as defined in (3.2), between a bound and an unbound quark-antiquark 
pair moving through the isotropic plasma (2.3) with velocities (from the rightmost curve to the 
leftmost curve) v = 0,0.35,0.85,0.996. The dipole is oriented orthogonally to its velocity. For 

V < Vtra,ns ouB has < imax, whcrcas for V > Wtrana ouc fiuds L = Lniaxj whcrc Wtrans — 0.45 Is thc 
transition velocity between the two behaviors. At v = the screening length and the maximum 
separation are ~ 0.24/r and Lmax — 0.27/r, respectively. 

3 Preliminaries 

In this paper we define the screening length Lg as the separation between a quark and 
an antiquark such that for i < Lg {i > Lg) it is energetically favorable for the quark- 
antiquark pair to be bound (unbound) [32, 33]. Obviously this satisfies Lg < L^^^, where 
L^ax is the maximum separation Lj^ax for which a bound quark-antiquark solution exists. 
We will determine Lg by comparing the action S{i) of the bound pair, which is a function 
of the quark-antiquark separation £, to the action ^unbound of the unbound system, i.e. by 
computing: 

AS{£) = S{e) - . (3.1) 

The screening length is the maximum value of i for which AS is positive (since we will 
work in Lorentzian signature). This may correspond to the value of i at which AS" crosses 
zero, in which case Lg < L^^^, or the maximum value of i for which a bound state exists, 
in which case Lg = L^^^. In the Euclidean version of our calculations, this criterion cor- 
responds to determining which configuration has the lowest free energy, which is therefore 
the configuration that is thermodynamically preferred. As shown in Fig. 4, for a meson 
moving through the isotropic plasma (2.3) one has Lg < L^-^^^ for v < f trans , whereas for 

V > Vtrans ouc fiuds that Lg = L„,ax5 where Wtrans — 0.45 is the transition velocity between 
the two behaviors [34-36] . These qualitative features extend to the anisotropic case, as we 
have illustrated in Fig. 5. The transition velocity decreases with the anisotropy, so for large 
a/T one has Lg = L^^x except for very low velocities. Similarly, if the ultra-relativistic 
limit v — )• 1 is taken at fixed a and T, then obviously v > Vtrans and again Lg = L^^^. 

All our calculations will be done in the rest frame of the quark-antiquark pair, to 
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Figure 5. Energy difference in an anisotropic plasma, as defined in (3.2), between a bound and 
an unbound quark-antiquark pair oriented along the transverse direction x and moving along the 
anisotropic direction z. All the curves on the left correspond to a/T = 12.2 and different velocities 
(from the rightmost curve to the leftmost curve) v = 0, 0.35, 0.85, 0.996. All the curves on the right 
correspond to the same velocity v = 0.25 and different anisotropics (from the rightmost curve to the 
leftmost curve) a/T — 0,6.5,43,744. For these anisotropies the corresponding transition velocities 
are respectively given by Wtrans — 0.45, 0.29, 0.19, 0.11. 



which we will refer as the dipole rest frame. Since any observable can be easily translated 
between this frame and the plasma rest frame, we will speak interchangeably of 'mesons in 
a plasma wind' and of 'mesons in motion in the plasma'. We emphasize however that all 
the physical quantities that we will present, e.g. the screening length, are computed in the 
dipole rest frame. 

The actions are scalar quantities, so A5dipoio = ASpia^ma- Moreover, in the dipole rest 
frame we have 

A'S'dipoio = "T'AE'dipol;, , (3.2) 

since the dipole is static in its own rest frame. In this expression i^dipoic is the energy 
(as opposed to the free energy) of the configuration and T = J dt is the length of the 
integration region in time. Thus we see that our criterion, which is based on comparing 
the actions, can also be thought of as a comparison between the energies of the bound and 
the unbound configurations in the dipole rest frame. 

We will see that the ultraviolet divergences in the string action associated to integrating 
all the way to the boundary of AdS cancel out in the difference (3.1), and neither the bound 
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nor the unbound actions possess infrared divergences associated to integrating all the way 
down to the horizon. This can be verified explicitly and it also follows from their relation 
to the energy in the rest frame of the dipole: While the energy of the unbound string 
pair possesses an infrared logarithmic divergence in the plasma rest frame [37], no such 
divergence is present in the dipole rest frame (see e.g. the discussion in [36]). 



4 Static dipole in an anisotropic plasma 

In an anisotropic plasma the screening length depends on the relative orientation between 
the dipole and the anisotropic direction z. Given the rotational symmetry in the xy-plane 
we assume without loss of generality that the dipole lies in the xz-plane, at an angle Q with 
the z-axis. We thus choose the static gauge t = T,a = u and specify the string embedding 
as 

X — )• sin x{u) , z — )■ cos 6 z{u) . (4-1) 
The string action takes the form 

t2 r /"Umax 1 I 

S = 72 dt du^J B (1 + Tncos"^ e z'^ + Tsm^ 9 x''^) , (4.2) 

2vra' J Jq u^^ 

where the 2 comes from the two branches of the string and ti^ax will be determined below. 
The conserved momenta associated to translation invariance in the x, z directions are given 

by 

1 dC BTsinOx' 

nx = — ^ = , (4.3) 

sm 8dx y2. Bh+FU cos2 6 z'^ + F sin^ 6 x'2) 



1 dC BFTicosOz' 

= = y — ■ (4.4) 

cos &dz u\ B(l+Fn cos2 ez'^+F sin^ 6 x'2) 



Inverting these relations we find 

, VTi CSC 9 u'^ Hx , scc6v?Iiz 

Vf^/bfu - (n2 + 71111)' ^fru^BTU - (n2 + u n2.) ' 

Substituting back in the action we arrive at 

^2 /■ /-Umax 1 B^TH 

2 dt du^ (4.6) 



2^a' J Jo ./BTH^^Hjl^+nW) 

For a U-shaped string describing a bound quark-antiquark pair the turning point u^^^ 
is determined in terms of the momenta by the condition that x'(u^^^) = z'{u^^^) — )• oo, i.e. 

by 

B:Fn-u'{ul + nul)l^_ = o. (4.7) 
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In turn, the momenta are determined by the boundary conditions that require the string 
endpoints to he a distance I apart from each other: 

dux' = \ du z' . (4.8) 



These two equations, together with (4.7), can be solved numerically to express the momenta 
and Uniax in terms of i. In this way the on-shell action (4.6) for a bound pair becomes a 
function of I alone. In order to determine Ls we subtract from this action the action of a 
static, unbound quark-antiquark pair, which is described by two straight strings hanging 
down from the boundary to the horizon. The action of this unbound pair is equal to (4.6) 
with the momenta set to zero and the range of integration extended down to the horizon: 



2,2 r run 

^unbound = 2 J dt J . (4.9) 

We obtain the screening length by numerically determining the value of i at which the 
difference S{£) — S'unbound crosses zero, since in the static case we always have Lg < ^max- 
The result for this difference as a function of i in the isotropic plasma [32, 33] described 
by eqn. (2.3) is plotted in Fig. 4, from which we see that the screening length is 

24 

Li,„(r) ~ — [static dipole] . (4.10) 

The scaling with the temperature is expected on dimensional grounds. In the isotropic case 
the temperature and the entropy density are related simply through (2.4), so this result 
can be recast as 

/ ^2^2 X 1/3 

Li,„(s) ~ 0.24 ( -^-^ j [static dipole] , (4.11) 

which will be useful later. 

The results in the anisotropic case are plotted in Figs. 6 and 7. Fig. 6 shows the 
screening length, for several orientations of the dipole, as a function of the anisotropy 
measured in units of the temperature (left) and the entropy density (right). The reason 
for working with both normalizations is that we wish to compare the screening length 
in the anisotropic plasma to that in the isotropic plasma, and this can be done at least 
in two different ways: the two plasmas can be taken to have the same temperatures but 
different entropy densities, or the same entropy densities but different temperatures. Fig. 7 
shows the screening length as a function of the dipole orientation for several values of the 
anisotropy. 

We see from Fig. 6(left) that Lg decreases monotonically as a increases, for any dipole 
orientation, if the temperature is kept fixed. We also see from Fig. 7(left) that this effect 
is more pronounced for a dipole oriented along the anisotropic direction. In contrast, 
the behavior of the screening length at constant entropy density depends on the dipole's 
orientation, as shown in Figs. 6(right) and 7(right). For dipole's aligned sufficiently close 
to the anisotropic direction the screening length decreases with the anisotropy, whereas for 
orientations sufficiently close to the transverse plane the screening length increases with 
the anisotropy. 
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Figure 6. Screening length as a function of the anisotropy for a static quark-antiquark dipole lying 
at an angle with the z-direction (from top to bottom) 9 = 7r/2, 7r/3, 7r/4, 7r/6, 0. The screening 
length is plotted in the appropriate units to facilitate comparison with the isotropic result for a 
plasma at the same temperature (left), or at the same entropy density (right). The isotropic result 
is given in eqs. (4.10) and (4.11). 
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Figure 7. Screening length for a quark-antiquark dipole lying at an angle 9 with the z-direction 
for anisotropics a/T — 12.2 (red, solid), 42.6 (maroon, coarsely dashed), 86 (violet, dashed), 744 
(orange, dot-dashed). The corresponding values in units of the entropy density are (in the same 
order) aN'^/^ / s^^'^ — 6.2,19,35,242. The screening length is plotted in the appropriate units to 
facilitate comparison with the isotropic result for a plasma at the same temperature (left), or at 
the same entropy density (right). The isotropic result is given in eqs. (4.10) and (4.11). 

5 Dipole in an anisotropic plasma wind 

In this section we will consider a static quark-antiquark pair in an anisotropic plasma that 
is moving with constant velocity with respect to the dipole — a dipole in an 'anisotropic 
plasma wind'. We will pay particular attention to the ultra-relativistic limit, which can be 
understood analytically.^ This limit, together with the static results from Sec. 4, will allow 
us to understand qualitatively the results at any velocity < w < 1. 

We will first rewrite the solution (2.2) in a boosted frame, and then place a dipole in 
^We recall that we first send the quark mass to infinity and then w — > 1 (see Sec. 1). 
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X 

Figure 8. Orientation of the dipole in an anisotropic plasma wind. The wind's velocity lies in the 
original xz-plane (before the boost (5.2)) at an angle 6y with respect to the anisotropic direction z. 
The quark lies at angles q = {x,y,z) = | (sin sin 1^9, sin cos 95, cos 0) with respect to the relabeled 
directions (after the boost (5.2)), and the antiquark lies at ~q. 

it — see Fig. 8. Given the rotational symmetry in the xy-plane we assume that the boost 
velocity is contained in the xz-plane, and that it lies at an angle 9^ with the z-axis. Thus 
we first rotate to a new coordinate system defined through 

t = i, 

X = z sin 9v + X cos 6^ , 

y = y, 

z = z cosOy — X sinOy , (5-1) 

and then perform a boost along the z-direction by setting 

i=j{t'-v z') , 
X = x' , 

y = y' , 

z = j[-vt' + z') , (5.2) 

where 7 = l/Vl — v'^ is the usual Lorentz factor. Below we will consider a dipole with an 
arbitrary orientation with respect to both the velocity of the plasma and the anisotropic 
direction z — see Fig. 8. We parametrize the orientation of the dipole by two angles 0, ip 
so that the quark hes at 

q = {x,y\z') = - {sin 9 simp, 9 cos if, cos 9) (5.3) 
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and the antiquark lies at —q. 

For notational simplicity, below we will drop the primes in the final set of coordinates. 
To avoid confusion, we emphasize that the direction 6^ of the plasma wind is always 
measured with respect to the original (x, y, z) axes, i.e. before the rotation and the boost 
above. In particular, motion within (outside) the transverse plane refers to a dipole in 
a plasma wind with 9^ = 7r/2 {6y ^ '?i"/2). In contrast, the orientation of the dipole 
is measured with respect to the final set of coordinates {x',y',z'). However, if instead 
of specifying the dipole's orientation through a pair (9, ip) we specify it by saying that 
the dipole is aligned with the x-, y- or z-directions, then we are referring to the original 
directions. Just as an illustration, consider the case of a plasma wind blowing along the 
original x-direction, i.e. a plasma wind with 9^ = 7r/2. Then we see from (5.1) and (5.2) 
that (x, z) ~ {z\ x'). Thus in this case by 'a dipole oriented along the x-direction' we mean 
a dipole with = 0. 

After dropping the primes from the final set of coordinates in (5.2) the five-dimensional 
part of the metric (2.2) takes the form 

ds^ = ^ [ -gttdt^ + gxxdx^ + dy^ + gzzdz"^ + gtxdt dx + gtzdt dz + gxzdx dz + —— 



where 



(5.4) 



BJ^ -v^isin^ 9^ +ncos'^ 9^) 

9tt = 2 ' 

1 — 

gxx = cos^ 9^+% sin^ 9^ , (5.6) 
sin^ 9^+% cos^ 9^ - v^BT 



oil! -r /LV^UD — U kJ J . . 

9zz = ^ ;2 ' ^^-^i 



gtx = i^i)jsin(20„), (5.8) 
V 1 — 



2v{BF -s\n^9^-'Hcos^9^) 

9tz = z — 2 ' y^-^) 

1 — f ^ 



gxz= ^^JL smi29y). (5.10) 

V 1 — 'V 

In order to determine the screening length for a generic velocity we need to compare 
the actions of a bound and an unbound quark-antiquark pair, as in the static case of 
Sec. 4. However, in the ultra-relativistic this is not strictly necessary because Lg = L^^^ 
(see Sec. 3). In other words, in this limit we only need to determine the maximum possible 
quark-antiquark separation for which a bound state exists. Nevertheless, for completeness 
we will briefly present the analysis of the unbound configuration. Each of the strings in the 
unbound pair is one of the trailing strings studied in [25], so the reader is referred to this 
reference for additional details. Note, however, that [25] worked in the plasma rest frame. 
Here we will work in the dipole's rest and focus on the ultra-relativistic limit. 
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5.1 Unbound quark-antiquark pair 



As in Sec. 4 we fix the static gauge t = t, a = u, and specify the embedding of the unbound 
string as 

X — )• x{u) , z — )• z{u) . (5-11) 

The embedding in the y-direction is simply y = because of rotational symmetry in the 
xy-plane and because the string is unbound. As we will see below, in the case of a bound 
string (dipole) the boundary conditions will generically imply a non-trivial embedding y{u). 



where 



The action for the unbound string reads 

L2 



L f /""^ 1 
— 2 J dt du-^^/F-^K^ + K,,x'2 + K,,z'^ + K^.x'z' , (5.12) 



K. 



Ko = gtt , 

BF{cos'^ + n sin2 6^) - Uv^ 



,2/1 , -7/„„„2, 



K,, = BT (sin^ 6^ + U cos^ 9^) , 
BT{1 - n) 



Introducing the conjugate momenta 

Br 

^''^^ unbound 



n. 



and solving for x' , z' we find 



dx' 



sin(2^tj 



9£„ 



where 



dz' 



1 



D = B-HTKo - u'' {K,,Ul + K^^Ul - K^zH^U,) 
Substituting into the action we arrive at 



27ra' 



u 



WD 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



The momenta are determined by the condition that (5.15) remain real for a string that 
extends all the way from the boundary to the horizon. Following [25] we analyze this 
condition by noting that D can be rewritten as 



27-4 

D = —N^N, - b 



(5.18) 
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where 



c= . A o • (5-19) 



— sin cos ' u 
As in [25] we must require that the zeros of the second summand in (5.18) coincide with 
one another and with those of and A^^^ . One of the zeros of the second summand occurs 
at a critical value u = Uc such that 

BcFc-dcv'^ = dc = Tic cos'^ 9^ + sm'^0v, (5.20) 

where Be = B{uc), etc. At this point we have 



AT ATI f^cosy„sm(/„ \^/Lc — i^jy^'j^^v^'-'-'-^v - ,^ 
^x^z\u. = — (^c - 1) dc dcUr, H , . (5.21 



{Tic — 1) cos 9v sin 0^ 



2 



Noting that He > 1 and that Kxz < 0, we see that D would be negative at Uc unless the 
momenta are related through 

^^= (^--^c)cos9^ ^ (5.22) 

dcVl^ 

Assuming this relation and requiring that the other zero in the second summand of (5.18) 
coincide with Uc yields 



2 _ BcJ'c dc 2 _ BcTcjy-c -if cos^ 9y sin^ 9y 

- , - ut{l-v^)dc ■ ^ ' 



Note that Hz does not vanish for any value of 9^, whereas Hx vanishes if 9^ = 0,7r/2. 
The reason is that for these two particular orientations the plasma wind blows along the 
original z- or x-directions and the string orients itself with the corresponding axis [25]. As 
a consequence, the momentum along the orthogonal axis vanishes. However, the changes 
of coordinates (5.1) and (5.2) always relabel the direction of motion as z, so after these 
changes the non-vanishing momentum is labelled 11^ irrespectively of whether 9y = oi 
9y = tt/2. 

We will analyze in detail the ultra- relativistic limit. This is facilitated by explicitly 
distinguishing the case of motion outside the transverse plane {9y ^ ''r/2) and motion 
within the transverse plane {9v = vr/2). 

5.1.1 Ultra-relativistic motion outside the transverse plane 

In the ultra-relativistic limit Uc approaches the boundary, i.e. Uc — ?■ 0, and we can use the 
near-boundary expansion (2.7) to determine it. The condition (5.20) yields in this limit 
[25] 

^ 1^ ~fj [Ov + ^/2] , (5.24) 

cos^ 9^ 

which when substituted in (5.23) gives the momenta 

2 g"^ cos^ 9y 2 cos^ 9^ sin^ 9^ 

^ 16(1 -^;2)2 ' ^ 16(1 -z;2) ' ^^'^^^ 
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In these expressions we have ignored subleading terms in an expansion in 1— f^, for example 
we have set ?; ~ 1, l-Lc ~ 1; etc. Note that in this expansion Hx is subleading with respect 

to n^. 

For later use we must evaluate how S'unbound scales with 1 — in the limit u — )• 1 . For 
this purpose we split the integration region, and hence the action (5.17), as 

C _ C<(1) _|_ 0(2) (r- r^a\ 

'--'unbound '-^unbound ^ '-^unbound ; 

where S^2,ovLnd is the action with the integral in u ranging between and Uc, and S'inLund is 
the action with the integral in u ranging between Uc and Uh- The reason for this separation 
is that in the first interval u is small and hence we will be able to use the near-boundary 
expressions (2.7), (5.24) and (5.25). In order to exhibit the dependence on 1 — of sf^bound 
explicitly, it is convenient to work with a rescaled variable r which remains finite in the 
V — >• 1 limit, defined though 



u 



-\/ 1 — v"^ ^ Uc = rc\/ 1 — . (5.27) 
In terms of this variable we get 

(1) 2 f p l-iaVcos20, + ... 

'-'unbound — ~ • ' I I Ur 



2Tia' Vl - 7 Jq ^2^1 _ ia2^2 cos2 - ^a^r^ cos^ 6^ + ... 

(5.28) 

The divergence near r = will cancel out with that in the action for the bound string. The 
integrand is smooth across r = re- The crucial point is that the result is O [(1 — f 2)"^/^] 
in the counting in powers of 1 — v^, and we will find this same scaling in the bound string 
action (see below). In contrast, 'S'^nLund scales as 1 — in the ultra-relativistic limit. The 
reason is that u is not small in units oi 1 — in the corresponding region of integration, 
so all the dependence comes from the fact that the action (5.17) scales as l/Hz ^ 1 — v"^ 
in this region. 

5.1.2 Ultra-relativistic motion vi^ithin the transverse plane 

In this case 9^ = 7r/2 and hence we see from (5.23) that Hx = 0. The condition (5.20) now 
gives [25] 

^c^y^^, (5.29) 

where 

121 

C = — - -F4 - ^4 , (5.30) 
57d 

and we recall that J-4, B4 are the coefficients that enter the near-boundary expansion (2.7). 
Substituting (5.29) into (5.23) and dropping subleading terms as before we obtain the 
momentum in the z-direction (recall that this corresponds to the original x-direction) : 



It is now convenient to work with a rescaled radial coordinate r defined through 

u = r{l-v^y/\ (5.32) 

SpHtting the unbound string action as before, we find 

q(i) - 2 f p + 

- W (1 - z;2)i/4 y rV l - 2Cr4 + . = " ^^"^^^ 

Again, the divergence near r = will cancel out with that in the action for the bound string, 
which will also be of O [(1 — ?;^)~^/^] in the counting in powers of 1 — u ^ (see below). In 
contrast, S'fnbound scales as l/Hx ~ — v'^ in the ultra-relativistic limit, and is therefore 
subleading. 

In summary, we find that in the ultra-relativistic limit 

O [(1 — w^)"^/^] if 7^ 7r/2 [outside the transverse plane] 



S., 



(5.34) 

O [(1 — w^)~^/^] if 9^ = tt/2 [within the transverse plane] . 



5.2 Bound quark-antiquark pair 

We now consider a dipole with an arbitrary orientation with respect to both the velocity 
of the plasma and the anisotropic direction z — see Fig. 8. As before we fix the static 
gauge T = t,a = u and specify the string embedding via three functions {x{u), y{u), z{u)) 
subject to the boundary conditions 

x'du , 







- sm a sm ip 
2 ^ 

- sin 6 cos (/9 = / y'du , 
2 Jo 

-COS 6 = / z'du, (5.35) 







where u^^^ is the turning point of the U-shaped string. The integral in the action of the 
bound string extends only up to this point and now includes a term proportional to y'^: 

r2 /• /■Umax 1 , 

S = - — 2 dt du^ ^J'-^Ko + K^^x'^ + Kyyv'^ + K,,z'^ + K^.x'z' . 

(5.36) 

All the K^s were defined in (5.13) except for Kyy, which is given by 

^yy - 1 -v'^ ' 

The momenta are defined as 

^ dC ^ dC ^ d£ , , 
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Inverting these equations we get 



y = > 



( \ \ 
z = , ^ — Kxz^x + Kxx^^ , (5.39) 

where 

D = BUFKq - {K,, lil + BFU + K^x - K^z Ho^n,) . (5.40) 
Substituting these expressions into the action (5.36) we get 

S = -^,2\,tr-,u^^. (5.41) 

As in the case of the unbound string, we will now distinguish between the cases of motion 
outside and within the transverse plane, focusing on the ultra-relativistic limit. 

5.2.1 Ultra-relativistic motion outside the transverse plane 

The turn-around point u^ax is defined by the condition D{u^^^) = 0. In the ultra-relativistic 
limit we expect that this point approaches the boundary for the string solution of interest, 
as in the isotropic case. Thus in this limit Umax can be determined by using the near- 
boundary expansions of the metric functions (2.7). 

In the limit u — )• we find the following expansions: 

K,,^l + + ■■■ , (5.42) 

a2«2sin6l^cos(9„ 
Kxz — , , (5.43) 

Kxx^l- — ^ + • • • , (5.44) 



4(1 -t;2) 

from which it follows that 

a^u^ cos^ 0. 



I?:^i-^^^^^f^-u^(n2+n2 + n2) + ... . (5.45) 



Similarly, the boundary conditions (5.35) take the form 



-sin0sin99~ / du—^Ilx + ---, (5.46) 
2 Jo V D 



- sin 9 cos — du —= liy + 

2 Jo VD 



I rUm.. ^2 / ^2^2 ^Qg2 
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In the ultra-relativistic limit, all the terms that we have omitted in the equations above, 
in particular in (5.45) and (5.47), are subleading with respect to the terms that we have 
retained provided the radial coordinate and the momenta scale as 

u = r^/\-v'^, Ui = -^, (5.47) 

1 — 

where r and pi are kept fixed in the limit f — t- 1. In terms of these rescaled variables (5.47) 
the boundary conditions (5.47) take the form 



- sm U sm If 



\/l - V^Pc,l2{p, 6v) , 



£ I 

- sin d cos (/? ~ \/l — v'^ Pyl2{p,dv 



2 

where the integral 



cose ^ Vl-v^Pz iMP,0.) ^Mp,0.)] , (5.48) 



Tn{p,0^)= / dr , (5.49) 



4 



cos^ Ov — r^{p1. +Py+ p1) 
is of 0(1) in the counting in powers in (1 — v"^), and is finite if n > 0. Further noting that 



a^u^cos^^^ . a?r'^ cos^Oy 



Ko = l-^^^^-^ + 0{u-)^l (5.50) 



S ^ -^,^=L=[z.2{pM-'^^^^Z,{p,e,)j j dt. (5.51) 



we see that the bound action scales as 

Since both this bound action and the unbound action (5.28) scale as (1 — u2)-i/2^ ^^j^g 
divergence at r = in the bound action coming from the I-2{p, 9^) integral would exactly 
cancel that in the unbound action in the difference (3.1). Moreover, by comparing the two 
actions we would conclude that the momenta pi introduced in (5.47) are indeed of 0(1) in 
the counting in powers of (1 — v'^) in the ultra-relativistic limit. It would then follow that 
the integrals Xn{p,Ov) are also of 0(1), and therefore that the screening length scales as 
~ (1 — w^)^/^ in the ultra-relativistic limit. However, as explained below (5.10), in the 
ultra-relativistic Lg = L^^^ is simply the maximum possible separation between a bound 
quark-antiquark pair, so it can be determined by maximizing i in (5.48) with respect to 
the momenta. Since the integrals are bounded from above for any value of the pi, and the 
maximum is u-independent, it follows that Ls = Lmax ~ (1 — t;^)!/^. 

5.2.2 Ultra-relativistic motion vi^ithin the transverse plane 

In this case 9^ = 7r/2 and the expansions of D and of the boundary conditions (5.35) 
become 

Dc^l-- - u\ul + Ul + nl) + --- (5.52) 

1 — y 



-18- 



and 



- sm sm 93 — / duu — -j^^^^^^^^=^^^^^^^^= + 

2 ^0 ^i-f^-n4(n2 + n2 + n2) 

- sm V cos if — duu — + 

2 /i_iz^j_^4(n2 + n2 + n2) 



cos ~ / — ^ ^ + • 

2 ^ ./i-ai-n4(n2 + n2 + n2) 



where C was defined in (5.30). As in the previous section, in the ultra-relativistic hmit all 
the terms that we have omitted in the equations above are subleading with respect to the 
terms that we have retained provided the radial coordinate and the momenta scale in this 
case as 

u = r{l-v')'/\ n, = -JL=, (5.53) 

V 1 — f 

where r and pt are kept fixed in the limit u — )• 1. In terms of the rescaled variables the 
boundary conditions (5.53) become 

^ sin 6* sin (1 - v"^)^/^ J2{p) , 

^sin6'cos(^~ {I - v'^y/'^PyJ2{p) , 

^cos^ ~ (1 -t;2)V4p^ (j2(p) - CMp)) , (5.54) 



where the integral 

Jn{p)= I dr (5.55) 
ll-r\C + pl+pl+pl) 

is of 0(1) in the counting in powers in (1 — w^), and is finite if n > 0. Further noting that 

Ko = l + Oiu^) ~ 1 - , (5.56) 

1 — t;^ 

we see that the bound action becomes 

' 27rQ' ( 



S^~T^^;^,[J-2{v)-CMp)) I dt. (5.57) 



Since both this bound action and the unbound action (5.33) scale as (1 — v^)~^/^, the 
divergence at r = in the bound action coming from the J'-2{p) integral would exactly 
cancel that in the unbound action in the difference (3.1). Moreover, by comparing the two 
actions we would conclude that the momenta pi introduced in (5.53) are indeed of 0(1) 
in the counting in powers of (1 — v^) in the ultra-relativistic limit. It would then follow 
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that the integrals Jn{p) are also of 0(1), and therefore that the screening length scales as 
~ (1 — f^)^/^ in the ultra-relativistic limit. However, as explained below (5.10), in the 
ultra-relativistic Lg = L^^^ is simply the maximum possible separation between a bound 
quark-antiquark pair, so it can be determined by maximizing i in (5.54) with respect to 
the momenta. Since the integrals are bounded from above for any value of the pi, and the 
maximum is ^-independent, it follows that Ls = Lmax ~ (1 — v"^)^/^. 

In summary, we conclude that in the dipole rest frame the screening length scales in 
the ultra-relativistic limit as 

{(1 — u^)-"^/^ if 6t, 7^ tt/2 [motion outside the transverse plane] 
(5.58) 
(1 — v'^)^^^ if 6^, = tt/2 [motion within the transverse plane] 

irrespectively of the dipole orientation. 

5.3 Isotropic limit 

The results above reduce to the isotropic result of Ref. [34, 35] in the limit a ^ 0. This 
limit is most easily recovered from the results for motion within the transverse plane, since 
some of the terms in the expansions in Section 5.2.1 vanish if a = 0, thus invalidating the 
analysis. In contrast, setting a = in Section 5.2.2 boils down to simply setting C to its 
isotropic value, which from (5.30) and (2.3) is 

C = -T4 = ^= vr^T^ . (5.59) 
K 

Since the value of C does not affect the ultra-relativistic scaling of the screening length, 
we recover the scaling 

Lisa ~ (1 — v'^)^^^ [isotropic plasma] (5.60) 

found in the isotropic case by the authors of [34, 35]. As in the anisotropic case, the 
ultra-relativistic scaling of the screening length is independent of the dipole's orientation. 
In fact, even for v < 1, the isotropic screening length depends only mildly on the dipole's 
orientation, as shown in Fig. 9. 

5.4 Numerical results for generic velocities 

Away from the ultra-relativistic limit the screening length must be obtained numerically. 
For this reason we have focused on a few representative cases, namely those in which both 
the direction of the plasma wind and the dipole's orientation are aligned with one of the 
original x, y, or z axes. Given the rotational symmetry in the xy- plane, there are only five 
inequivalent cases to consider, because if the wind 'blows' in the z-direction then orienting 
the dipole along x or y gives identical physics. In each case, we plot the screening length 
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Figure 9. Screening length for a dipole moving through an isotropic plasma in a direction orthog- 
onal (top, blue curve) or parallel (bottom, orange curve) to its orientation. 

both as a function of the velocity v for different degrees of anisotropy a, and also as a 
function of the degree of anisotropy for different values of the velocity. In each case the 
result can be qualitatively understood combining the static results from Sec. 4 and the 
ultra-relativistic behavior derived analytically in Section 5. We recall that in all cases 
below, by 'a dipole oriented along x, y or z' we are referring to the original directions 
before the rotation (5.1) and the boost (5.2). 

Wind along z and dipole along z. The numerical results are shown in Figs. 10 and 11. 
The curves in Fig. 10 start at v = with the same value as the 9 = static result shown in 
Fig. 7, and that they vanish as (1 — t;^)^/^ in the limit f — >• 1, in agreement with (5.58) (top 
line) and (5.60). The screening length decreases with the anisotropy, irrespectively of 
whether T or s are kept fixed. 

Wind along z and dipole along x. The numerical results are shown in Figs. 12 and 
13. We see that the curves in Fig. 12 start at v = with the same value as the 9 = 7r/2 
static result shown in Fig. 7, and that they vanish as (1 — 

^2)1/4 

the limit u — >• 1, in 

agreement with (5.58) (top line) and (5.60). In this case the screening length decreases with 
the anisotropy for any velocity provided the temperature is kept fixed. The same behavior 
is found at constant entropy density for high enough velocities, whereas for low velocities 
the screening length at constant s actually increases with a. 

Wind along x and dipole along x. The numerical results are shown in Figs. 14 and 
15. The curves in Fig. 14 start at f = with the same value as the 9 = Tr/2 static result 
shown in Fig. 7, and that they approach a finite, non-zero value as v — >• 1, in agreement 
with (5.58)(bottom line) and (5.60). As in previous cases, the screening length decreases 
with the anisotropy for any velocity provided the temperature is kept fixed. The opposite 
behavior is found at constant s. 

Wind along x and dipole along y. The numerical results are shown in Figs. 16 and 
17. We see that the curves in Fig. 16 start at v = with the same value as the 9 = 7r/2 
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Figure 10. Screening length for a plasma wind along the z-direction and a dipole oriented 
along the ^-direction, for four different values of the anisotropy (from top to bottom) a/T = 
12.2, 42.6, 86, 744. The corresponding values in units of the entropy density are (in the same 
order) aN'^/^/s^/'^ = 6.2,19,35,242. The screening length is plotted in the appropriate units to 
facilitate comparison with the isotropic result for a plasma at the same temperature (left), or at 
the same entropy density (right). The isotropic result is plotted in Fig. 9, and its ultra-relativistic 
behavior is given in eq. (5.60). At u = the curves agree with the 9 = Q values in Fig. 7. As w — 1 
they vanish as (1 — u^)^/^, in agreement with (5.58)(top line) and (5.60). 
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Figure 11. Screening length for a plasma wind along the z-direction and a dipole oriented along 
the z-direction, at five different velocities (from top to bottom) v — 0.25, 0.5, 0.7, 0.9, 0.9995. The 
screening length is plotted in the appropriate units to facilitate comparison with the isotropic result 
for a plasma at the same temperature (left), or at the same entropy density (right). The isotropic 
result is plotted in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). 



static result shown in Fig. 7, and that they approach a finite, non-zero value as u — )• 1, in 
agreement with (5.58) (bottom line) and (5.60). The qualitative behavior in as in the case 
of motion and orientation along x. 

Wind along x and dipole along z. The numerical results are shown in Figs. 18 and 
19. We see that the curves in Fig. 18 start at v = with the same value as the = 
static result shown in Fig. 7, and that they approach a finite, non-zero value as u — )• 1, in 
agreement with (5.58) (bottom line) and (5.60). The screening length decreases with the 
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Figure 12. Screening length for a plasma wind along the z-direction and a dipole oriented along 
the x-direction, for four different values of the anisotropy a/T = 12.2 (red, solid), 42.6 (maroon, 
coarsely dashed), 86 (violet, dashed), 744 (orange, dot-dashed). The corresponding values in units 
of the entropy density are (in the same order) aN"^/^ / s'^/'^ = 6.2, 19,35,242. The screening length 
is plotted in the appropriate units to facilitate comparison with the isotropic result for a plasma at 
the same temperature (left), or at the same entropy density (right). The isotropic result is plotted 
in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). At i; = the curves agree with 
the 6 = 7r/2 values in Fig. 7. As w — 1 they vanish as (1 — v^)^/^, in agreement with (5.58)(top 
line) and (5.60). 




a/T a/s^/"^ 

Figure 13. Screening length for a plasma wind along the z-direction and a dipole oriented along 
the x-direction, at five different velocities (from top to bottom) v = 0.25, 0.5, 0.7, 0.9, 0.9995. The 
screening length is plotted in the appropriate units to facilitate comparison with the isotropic result 
for a plasma at the same temperature (left), or at the same entropy density (right). The isotropic 
result is plotted in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). 

anisotropy for any velocity provided the temperature is kept fixed. The same is true at 
large anisotropics if the entropy density is kept fixed. 

6 Dissociation temperature and dissociation anisotropy 

In previous sections we have focused on computing the screening length in an anisotropic 
plasma, Ls{T,a), and on comparing it to its isotropic counterpart Li^„ = Ls(T,0). The 
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Figure 14. Screening length for a plasma wind along the a;-direction and a dipole oriented 
along the a:-direction, for four different values of the anisotropy (from top to bottom) a/T = 
12.2, 42.6, 86, 744. The corresponding values in units of the entropy density are (in the same 
order) aN"^/^ / s^^'^ = 6.2,19,35,242. The screening length is plotted in the appropriate units to 
facilitate comparison with the isotropic result for a plasma at the same temperature (left), or at 
the same entropy density (right). The isotropic result is plotted in Fig. 9, and its ultra-relativistic 
behavior is given in eq. (5.60). At w = the curves agree with the Q = 7r/2 values in Fig. 7. As 
w — >■ 1 they approach a finite, non-zero value, in agreement with (5.58) (bottom line) and (5.60). 
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Figure 15. Screening length for a plasma wind along the a;-direction and a dipole oriented along 
the x-direction, at five different velocities v =0.25 (yellow, dot-dashed), 0.5 (green, short dashed), 
0.7 (brown, medium dashed), 0.9 (cyan, long dashed), 0.9995 (blue, sohd). The screening length is 
plotted in the appropriate units to facilitate comparison with the isotropic result for a plasma at 
the same temperature (left), or at the same entropy density (right). The isotropic result is plotted 
in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). 



screening length characterizes the dissociation of a quark-antiquark pair for fixed T and 
a: A pair separated a distance i < Ls forms a bound state, but if i is increased above 
Lg then the bound state dissociates. Similarly, one may define a dissociation temperature 
TdissicL,^) that characterizes the dissociation of a quark-antiquark pair of fixed size £ in a 
plasma with a given degree of anisotropy a: for T < T^i^^ the pair forms a bound state, 
but if T is increased above T^i^s then the bound state dissociates. Analogously, one may 
define a dissociation anisotropy aiiss{T,£) such that a bound state forms for a < Odiss but 
not for a > a^i^a- It is useful to think of the three-dimensional space parametrized by 
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Figure 16. Screening length for a plasma wind along the a;-direction and a dipole oriented along 
the y-direction, for four different values of the anisotropy a/T = 12.2 (red, solid), 42.6 (maroon, 
coarsely dashed), 86 (violet, dashed), 744 (orange, dot-dashed). The corresponding values in units 
of the entropy density are (in the same order) aN^^^ /s-^^^ = 6.2, 19,35,242. The screening length 
is plotted in the appropriate units to facilitate comparison with the isotropic result for a plasma at 
the same temperature (left), or at the same entropy density (right). The isotropic result is plotted 
in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). At v = the curves agree with 
the 9 = IT 1 2 values in Fig. 7. As w — >■ 1 they approach a finite, non-zero value, in agreement with 
(5.58) (bottom line) and (5.60). 
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Figure 17. Screening length for a plasma wind along the a;-direction and a dipole oriented along 
the 2/-direction, at five difi'erent velocities v =0.25 (yellow, dot-dashed), 0.5 (green, short dashed), 
0.7 (brown, medium dashed), 0.9 (cyan, long dashed), 0.9995 (blue, sohd). The screening length is 
plotted in the appropriate units to facilitate comparison with the isotropic result for a plasma at 
the same temperature (left), or at the same entropy density (right). The isotropic result is plotted 
in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). 



(T, a,€) as divided in two disconnected regions by a two-dimensional surface: in one region 
quark-antiquark pairs bind together, while in the other one they do not. The functions 
Ls{T, a), Tdiss(a, £) and a^i^siT, i) are then simply different parametrizations of the dividing 
surface. It is therefore clear that if a triplet (T, a, I) lies on the dividing surface then 

TLs{a,T) = T^,,,{a,£)£, aLs{T,a) = a^UT,£)£ , etc. (6.1) 
In this section we will focus on the qualitative form of T^i^s and Odiss- As we will 
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Figure 18. Screening length for a plasma wind along the x-direction and a dipole oriented 
along the ^-direction, for four different values of the anisotropy (from top to bottom) a/T — 
12.2, 42.6, 86, 744. The corresponding values in units of the entropy density are (in the same 
order) aN'^/^/s^/'^ = 6.2,19,35,242. The screening length is plotted in the appropriate units to 
facilitate comparison with the isotropic result for a plasma at the same temperature (left), or at 
the same entropy density (right). The isotropic result is plotted in Fig. 9, and its ultra-relativistic 
behavior is given in eq. (5.60). At v — the curves agree with the 6 = Q values in Fig. 7. As w — 1 
they approach a finite, non-zero value, in agreement with (5.58) (bottom line) and (5.60). 
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Figure 19. Screening length for a plasma wind along the x-direction and a dipole oriented along 
the z-direction, at five different velocities (from bottom to top) v ~ 0.25, 0.5, 0.7, 0.9, 0.9995. The 
screening length is plotted in the appropriate units to facilitate comparison with the isotropic result 
for a plasma at the same temperature (left), or at the same entropy density (right). The isotropic 
result is plotted in Fig. 9, and its ultra-relativistic behavior is given in eq. (5.60). 



see, most of the analysis follows from the asymptotic behavior of the screening length for 
a ^ T. This means that, at the qualitative level, most of the results that we will obtain 
would also apply if we were to replace the temperature by the entropy density as one of our 
variables. The reason is that, by virtue of (2.5), the limit a ^ T corresponds to the limit 
a ^ s^^^ and vice versa. In addition, we will see that for generic dipole's orientations and 
velocities, the large- anisotropy limit is entirely controlled by the near-boundary behavior 
of the metric at 0{u^), which depends solely on a and is therefore completely insensitive 
to the values of the temperature or of the entropy density. 
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Figure 20. Proper velocity in the z-direction at a position u away from the boundary, as defined 
in (6.2), for different values of a/T. From right to left, a/T = f .38, 33, 86, 249. 

The key point in the large-a analysis is the requirement that no point on the string 
can move faster than the local speed of light in the bulk. Consider a meson moving with a 
velocity v that has a non-zero component Vz along the z-direction. Then we see from (2.2) 
that the proper velocity along this direction of a point on the string sitting at a value u of 
the radial coordinate is 

VpropcAU) = Vzi —- = Vzi . . , . . (6.2) 

V 9tt{u) y F[u)B[u) 

The function T-L{u) increases monotonically from the boundary to the horizon, and is does 
so more steeply as a/T increases, as illustrated in Fig. 1. The combination T{u)B{u) has 
the opposite behavior, as expected from the fact that gravity is attractive: it decreases 
monotonically from the boundary to the horizon. In the isotropic case T-L = 1 and J-B 
decreases more steeply as T increases. This is thus the first hint that increasing the 
anisotropy has an effect similar to increasing the temperature: both make f proper (^i) a more 
steeply increasing function away from the boundary. We have illustrated the effect of the 
anisotropy in Fig. 20, where we see that Wpropcr/^z becomes a steeper function of u as a/T 
increases. 

It follows that, for fixed Vz ^ 0, there is a maximum value of timax beyond which 
^'proper becomes superluminal, so no string solution can penetrate to u > Umax- As we 
will corroborate numerically, this upper bound on u^^^ translates into an upper bound 
on Lg. Moreover, Wmax decreases as a/T increases. This means that for sufficiently large 
anisotropics we can use the near-boundary expansions (2.7) in order to determine Lg, in 
analogy to what we did in the ultra-relativistic limit. As in that case, for Vz ^ the analysis 
is controlled by the 0{u'^) terms in (2.7). The key point is that these terms depend on a 
but not on T, so by dimensional analysis it follows that timax ~ cl~^ and Ls ~ a~^ in the 
limit a/T ^ 1. This limit can be understood as a — t- oo at fixed T, or as T — t- at fixed 
a. We thus conclude that, even at T = 0, a generic meson will dissociate for a sufficiently 
large anisotropy a^i^s- 
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Mesons at rest and mesons whose velocity is exactly aligned with the transverse plane 
constitute an exception to the argument above, since in this case Vz = and their physics 
is mostly insensitive to the function T-L{u) which characterizes the anisotropic direction. 
Therefore in this case we expect that Umax and Lg will remain finite as we send a — >■ oo at 
fixed T, and hence that dimensional analysis will imply Lg ~ T^^. 

In summary, the heuristic argument above suggests that in the limit a/T ^ 1 we 
should have 



Ls{T,a) ~ <^ 



const. X T~ if the meson is static or in motion within the transverse plane, 
const. X a^^ otherwise. 



(6.3) 



The constants may depend on all the dimensionless parameters such as the velocity and the 
dipole's orientation. We will refer to the behavior in the second line as 'generic' and to that 
in the first line as 'non-generic', since the latter only applies if the velocity is exactly zero or 
if the motion is exactly aligned with the transverse plane. The generic behavior is of course 
consistent with the analysis of Sec. 5.2.1. Indeed, we saw in that section that for motion 
outside the transverse plane the ultra-relativistic behavior of Lg is entirely controlled by 
the O(ti^) terms in the metric, which depend on a but not on T. 

Fig. 21 shows our numerical results for u^^^, in units of and a^^, as a function 
of a/T, for the five physically distinct cases discussed in Sec. 5.4. From the continuous, 
magenta curves in the first two rows we see that n^ax goes to zero at large a/T in the cases 
of motion along z, irrespectively of the dipole's orientation. In contrast, we see that li^ax 
does not go zero for a static meson (dashed, blue curves) or for a meson moving along the 
x-direction (continuous, magenta curves in the last three rows). 

Recalling that the isotropic screening length is of the form Li^^ oc 1/T, we see that 
the quantity plotted on the vertical axes in Figs. 6, 11, 13, 15, 17 and 19 is precisely 
proportional to TLs{T,a). However, the asymptotic behavior (6.3) is not apparent in 
these plots because in most cases the horizontal axes do not extend to high enough values 
of a/T. For this reason we have illustrated the two possible asymptotic behaviors of Lg in 
Fig. 22, where we have extended the horizontal axes to larger values of a/T. We see from 
the continuous, magenta curves in the first two rows that ~ 1/a for motion along the 
z-direction. For motion within the transverse plane we see from the same curves in the last 
three rows that Lg ~ 1/T. This approximate scaling relation seems to hold quite precisely 
for a dipole oriented within the transverse plane (3rd and 4th rows), whereas for a dipole 
oriented in the z-direction the product TLg seems to retain a slight (perhaps logarithmic) 
dependence on a/T at large a/T. We can draw similar conclusions from the dashed, blue 
curves in the figure, which correspond to static mesons. We see that for mesons oriented 
within the transverse plane (2nd, 3rd and 4th rows) the relation TLg ~ constant holds 
quite precisely, whereas for mesons oriented in the z-direction (1st and 5th rows) there 
seems to be some slight residual dependence on a/T at large a/T. 
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Figure 21. Value of Umax in units of 1/T (left) or 1/a (right), as a function of the ratio a/T, for 
a dipole at rest (dashed, blue curve) and for a dipole moving with v = 0.45 (continuous, magenta 
curve). The first letter on the top right corner of each plot indicates the direction of motion, and 
the second one indicates the orientation of the dipole. 
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Figure 22. Screening length in units of 1/T (left) or 1/a (right), as a function of the ratio a/T, for 
a dipole at rest (dashed, blue curve) and for a dipole moving with v = 0.45 (continuous, magenta 
curve). The first letter on the top right corner of each plot indicates the direction of motion, and 
the second one indicates the orientation of the dipole. 
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Figure 23. Dissociation temperature (left) T^^^^{a^l) = i~^f{al) and dissociation anisotropy 
(right) adiss(r, ^) — £~^g{T£) for a dipole at rest (dashed curves) and for a dipolc moving along the 
z-direction with v = 0.45 (continuous curves). The orientation of the dipole is indicated by a letter 
next to each curve. 




Combining the two plots on the left and the right columns of Fig. 22 we can eliminate 
a/T and obtain TLs as a function of aLg and vice versa. Recalling (6.1) we see that we can 
interpret the result in the first case as T^iss{a,i) = £~^f{a£), whereas in the second case we 
get a^isa{T,i) = i~^g{Ti). The functions / and g are the curves shown in Fig. 23(left) and 
Fig. 23 (right), respectively. The right plot is of course the mirror image along a 45 degree 
line of the left plot. We see in Fig. 23(left) that the dissociation temperature decreases 
monotonically with increasing anisotropy and vanishes at a£ ~ 9.75 (for the chosen velocity 
and orientation). On the right plot this corresponds to the dissociation anisotropy at zero 
temperature. As anticipated above, even at zero temperature, a generic meson of size i will 
dissociate if the anisotropy is increased above anssiT = 0,i) oc The proportionality 
constant in this relation is a decreasing function of the meson velocity in the plasma. This 
implies that for a fixed anisotropy there is a limiting velocity Vn^^ above which a meson will 
dissociate, even at zero temperature. The form of Vii^{ai) for T = is plotted in Fig. 24. 

The existence of a limiting velocity for quarkonium mesons is well known in a strongly 
coupled isotropic plasma [38, 39], in which case the dissociation at f = is caused by 
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Figure 25. Dissociation temperature for a meson moving along the x-direction and oriented along 
the z-direction (left) or along the x-direction (right). Each curve corresponds to a fixed value of 
the product a£ = (blue curve), 1.4 (green curve), 25 (red curve). 



the temperature. What we see here is that in our anisotropic plasma this behavior persists 
as T — )• for generic motion. In this limit it is the anisotropy that is responsible for the 
dissociation. In the case of ultra-relativistic motion the relation between a^i^s or T^i^s and 
Vii^ can be obtained by combining the scalings (5.58) and (6.3). For generic motion these 
relations yield 

a,UT, ^) ~ ^(1 - viy/^ , [a » T , v,,^, < 1] (6.4) 

whereas for motion within the transverse plane we obtain 

T,Ua, ^) ~ ^(1 - viy/* . [a > T , v,,,^ < 1] (6.5) 

The scaling (6.5) agrees with the isotropic result [34, 35] and illustrates the fact that, 
for motion within the transverse plane, the limiting velocity in our anisotropic plasma 
approaches unity as T — )• 0. This behavior is the same for a meson at rest, as illustrated in 
Fig. 23, where we see that a sufficiently small meson will remain bound in the plasma for 
any value of the anisotropy provided the plasma is cold enough. In fact, the form of the 
dissociation temperature for all anisotropies and all velocities within the transverse plane 
is qualitatively analogous to that of the isotropic shown in Fig. 25. The fact that 

the curves in this figure approximately overlap one another signals that the dependence of 
the dissociation temperature on v and ai can be approximately factorized over the entire 
range < v < 1. 

In contrast, for generic motion we saw above that the limiting velocity is subluminal 
even at T = 0, Viin^{T = 0,a£) < 1. Increasing the temperature simply decreases the value 
of the limiting velocity, Vu^iTl, a£) < Vii^{T = 0, a£). Turning these statements around we 
see that, at a fixed anisotropy, the dissociation temperature is a decreasing function of 
the velocity that vanishes at v = Viim(^ = 0,a£). This is illustrated in Fig. 26, where we 
see that Vii^{T = Q,a£) decreases as the anisotropy increases, in agreement with Fig. 24. 
In order to facilitate comparison with the isotropic results of [34-36], in Fig. 26 we have 
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Figure 26. Dissociation temperature for a meson moving along the z-direction and oriented along 
the a;-direction (left) or along the z-direction (right). Each curve corresponds to a fixed value of 
the product a^. From right to left, a£ = 0, 1, 5.4, 25. 



chosen to normalize the dissociation temperature by its value at f = instead of by the 
dipole's size Our numerical results suggest that as v approaches Vi^^ the dissociation 
temperature may vanish as 

In this equation v^^^ = Vii^{T = 0,a£) and e = e{a£) > is an anisotropy-dependent 
exponent. Unfortunately, the limit v — )• v^^ is difficult to analyze numerically, so our 
results are not precise enough to allow us to establish (6.6) unambiguously. To emphasize 
this point, in Fig. 26 we have plotted as discontinuous the part of the curves between the 
last two data points. The last point lies on the horizontal axis at {v, T) = {vn^, 0), and the 
penultimate point lies at a certain height at {v < Vu^^jT > 0). Since this last bit of the 
curves is an interpolation between these data points, it is difficult to establish whether the 
slopes of the curves diverge as they meet the horizontal axis, as would be implied by the 
scaling (6.6). Presumably, this scaling could be verified or falsified analytically by including 
the first correction in T/a to the scaling in the second line of (6.3). 



7 Discussion 



We have considered an anisotropic = 4 SYM plasma in which the x, y directions are 
rotationally symmetric, but the z-direction is not. In the context of heavy ion collisions the 
latter would correspond to the beam direction, and the former to the transverse plane. The 
screening length of a quarkonium meson in motion in the plasma depends on the relative 
orientation between these directions, on the one hand, and the direction of motion of the 
meson and its orientation, on the other. This dependence can be parametrized by three 
angles {9^, 9, ip), as shown in Fig. 8. We have determined the screening length for the most 
general geometric parameters and for any anisotropy. Our results are valid in the strong- 
coupling, large- A^c limit, since we have obtained them by means of the gravity dual [19, 20] 
of the anisotropic = 4 plasma. The anisotropy is induced by a position-dependent theta 
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term in the gauge theory, or equivalently by a position-dependent axion on the gravity side. 
One may therefore wonder how sensitive the conclusions may be to the specific source of 
the anisotropy. In this respect it is useful to note that the gravity calculation involves 
only the coupling of the string to the background metric. This means that any anisotropy 
that gives rise to a qualitatively similar metric (and no Neveu-Schwarz i?-field) will yield 
qualitatively similar results for the screening length, irrespectively of the form of the rest 
of supergravity fields. 

An example of a rather robust conclusion is the ultra-relativistic behavior^ of the 
screening length (5.58), which for motion not exactly aligned with the transverse plane is 
Ls ~ (1 — v^Y^'^. The 1/2 exponent contrasts with the 1/4 isotropic result [34, 35], and 
follows from the fact that the near-boundary fall-off of the metric (2.2) takes the schematic 
form 



As V grows closer and closer to 1 the point of maximum penetration of the string into the 
bulk, Umax, moves closer and closer to the AdS boundary at n = 0. As a consequence, the 
physics in this limit is solely controlled by the near-boundary behavior of the metric. For 
generic motion the behavior is in fact governed by the O(u^) terms alone, and a simple 
scaling argument then leads to the 1/2 exponent above. In the isotropic case the O(u^) 
terms are absent and the same scaling argument leads to the 1/4 exponent. 

In fact, a similar reasoning allowed us to determine the large- anisotropy limit. Since 
the metric component gzz oc T-L{u) grows as one moves from the boundary to the horizon, 
a subluminal velocity of the meson at the boundary would eventually translate into a 
superluminal proper velocity (6.2) at a sufficiently large value of u.'^ This sets an upper 
limit on the maximum penetration length Umax of the string into the bulk and hence on Lg . 
Moreover, gzz becomes steeper as a/T increases, so in the limit a/T S> 1 the point u^ax 
approaches the AdS boundary (unless the motion is aligned with the transverse plane), just 
as in the ultra-relativistic limit. In this limit the physics is again controlled by the O(m^) 
terms in the metric, which depend on a but not on T. Therefore dimensional analysis 
implies that Lg = const, x a^^, were the proportionality 'constant' is a decreasing function 
of the velocity. This led us to one of our main conclusions: even in the limit T — )• 0, a 
generic meson of size £ will dissociate at some high enough anisotropy a^i^^ ^ . Similarly, 
for fixed o and T, even if T = 0, a generic meson will dissociate if its velocity exceeds a 
limiting velocity ^^(a, T) < 1, as shown in Fig. 24 for T = 0. As explained in Sec. 6, the 
conclusions in this paragraph would remain unchanged if we worked at constant entropy 
density instead of at constant temperature, since in the limit a ^ s^l^ the physics would 
again be controlled only by the O(u^) terms in the metric. 

The above discussion makes it clear that, at the qualitative level, much of the physics 
depends only on a few features of the solution: The presence of the gl^i, term in the near- 

^We recall that we first send the quark mass to infinity and then t; — > 1 (see Sec. 1). 
^Note that the overall conformal factor in (2.2) plays no role in this argument, since it cancels out 
in the ratio (6.2). 




(7.1) 
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boundary expansion of the metric, the fact that the metric (7.1) be non-boost-invariant 
at order u"^ (i.e. that gjjj not be proportional to r]^^), and the fact that gzz increases as 
a function of both u and a/T.^ The second condition is necessary because otherwise the 
physics of a meson in motion would be equivalent to that of a meson at rest, and we 
have seen that the latter is very similar to that of a meson in an isotropic plasma. The 
third condition ensures that Umax moves close to the boundary as a/T increases. Note that 
adding temperature to an otherwise boost-invariant metric will only affect g^fj , and thus 

(2) 

this is not enough to make guj non-boost-invariant. This conclusion is consistent with the 

(2) 

fact that gliJ is only a function of the external sources which the theory is coupled to. 

From the gauge theory viewpoint, some heuristic intuition can be gained by recalling 
that the anisotropy is induced by dissolving along the z-direction objects that extend 
along the xy-directions [19, 20, 29]. The number density of such objects along the z- 
direction, dn/dz, is proportional to a. On the gravity side these are D7-branes that wrap 
the five-sphere in the metric (2.2), extend along the xy-directions, and are homogeneously 
distributed in the z-direction. Increasing a has a large effect on the entropy density per 
unit 3- volume in the xyz-directions, in the sense that s/T^ — )• oo as a/T ^ oo, as shown 
in Fig. 2(left). In contrast, the entropy density per unit 2-area in the xy-directions on a 
constant-z slice, s^°/r^, approaches a constant in the limit a/T — )• oo. This is illustrated 
in Fig. 2(right), which is based on our numerical calculations, but it can also be proven 
analytically following the argument in Sec. 2.5 of Ref. [29]. In view of these differences, it is 
perhaps not surprising that the anisotropy has the largest effect on the physics of mesons 
moving along the z-direction, and the smallest effect on the physics of mesons moving 
within the transverse plane. Mesons at rest are also more sensitive to the anisotropy if 
they extend along the z-direction than if they are contained within the transverse plane. 
Presumably, the correct intuition behind this physics is that moving against the D7-branes 
is harder than moving along them. 

We close with a few comments on existing weak-coupling results on the physics of 
quarkonium dissociation in the real-world QGP in the presence of anisotropies. To the 
best of our knowledge no results at non-zero velocity exist, so we will limit ourselves to 
the static case. We emphasize though that any comparison between these results and ours 
should be interpreted with caution, because the sources of anisotropy in the QGP created 
in a heavy ion collision and in our system are different. In the QGP the anisotropy is 
dynamical in the sense that it is due to the initial distribution of particles in momentum 
space, which will evolve in time and eventually become isotropic. In contrast, in our case the 
anisotropy is due to an external source that keeps the system in an equilibrium anisotropic 
state that will not evolve in time. We hope that, nevertheless, our system might provide 
a good toy model for processes whose characteristic time scale is sufficiently shorter than 
the time scale controlling the time evolution of the QGP. 

A general conclusion of Refs. [40-42] is that, if the comparison between the anisotropic 
plasma and its isotropic counterpart is made at equal temperatures, then the screening 

*Again, up to possible overall conformal factors. 
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length increases with the anisotropy. This effect occurs for dipoles oriented both along and 
orthogonally to the anisotropic direction, but it is more pronounced for dipoles along the 
anisotropic direction. The dependence on the anisotropy in these weak-coupling results is 
the opposite of what we find in our strongly coupled plasma. In our case the screening 
length in the anisotropic plasma is smaller than in its isotropic counterpart if both plasmas 
are taken to have the same temperature, as shown in Fig. 6(left). We also find that the effect 
is more pronounced for dipoles extending along the anisotropic direction, as illustrated in 
Fig. 7(left). 

Refs. [42, 43] argued that if the comparison between the anisotropic and the isotropic 
plasmas is made at equal entropy densities, then the physics of quarkonium dissociation 
exhibits little or no sensitivity to the value of the anisotropy. This is again in contrast to 
our results since, as shown in Fig. 6(right) and Fig. 7(right), the screening length in this 
case is just as sensitive to the anisotropy as in the equal-temperature comparison. The 
difference in the equal-entropy case is simply that the screening length may increase or 
decrease with the anisotropy depending on the dipole's orientation. 
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